import numpy as np

def modindex(natural, noise):
	a = np.load(natural)
	b = np.load(noise)

	c = (a - b) / (a + b)

	return np.mean(c[np.where(np.isnan(c) == False)])
